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We calculate next-to-leading-order (NLO) corrections to exclusive processes in kr factorization 
theorem, taking 717* — » 7 as an example. Partons off-shell by k% are considered in both the 
quark diagrams from full QCD and the effective diagrams for the pion wave function. The gauge 
dependences in the above two sets of diagrams cancel, when deriving the fc-r-dependent hard kernel as 
their difference. The gauge invariance of the hard kernel is then proven to all orders by induction. 
The light-cone singularities in the fey-dependent pion wave function are regularized by rotating 
the Wilson lines away from the light cone. This regularization introduces a factorization-scheme 
dependence into the hard kernel, which can be minimized in the standard way. Both the large double 
logarithms In 2 kr and In 2 x, x being a parton momentum fraction, arise from the loop correction 
(3JT)' to the virtual photon vertex, the former being absorbed into the pion wave function and organized 

[""t , by the kr resummation, and the latter absorbed into a jet function and organized by the threshold 

resummation. The NLO corrections are found to be only few-percent for n-y* — > 7, if setting the 
factorization scale to the momentum transfer from the virtual photon. 
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p I. INTRODUCTION 

<D ' 

kx factorization theorem [3, [3, H, 0, H, @ , as a fundamental tool of perturbative QCD (PQCD), has been widely 
applied to inclusive and exclusive processes. It has been pointed out that kx factorization theorem is appropriate for 
processes dominated by contributions from small parton momentum factions x . Its application to exclusive B meson 
decays has led to the PQCD approach [1, [9|, Hoi. Illl. H2I ], which is free of the sing ularities from the end-point regions 
of x that usually appear in collinear factorization theorem [l3|, [3, EH [H, EE Il8| | Several aspects of kr factorization 
theorem have been studied. For example, a naive definition of /c^-dependent hadron wave functions, in which the 
(f-) ■ coordinate of a quark field is simply shifted by a transverse distance, contains light-cone divergences Modified 
definitions to remove these divergences have been proposed in [ltj 0, [2l| . The B meson wave function defined in 
kx factorization theorem is normalizable [2fj| . while the B meson distribution amplitude in collinear factorization 
theorem is not [HI, HH], when evolution effects are taken into account. The Sudakov resummation [il. IR l24l . I25I of the 
large double logarithm In 2 kr is essential for improving perturbative expansion in kx factorization theorem |26t HtJ ■ 
The current application of kx factorization theorem to exclusive processes is mainly made at leading order (LO) 
in the strong coupling constant a s [2q |: the important logarithms in hadron wave functions have been organized 
to all orders, but hard kernels are still evaluated at tree level. To demonstrate that kx factorization theorem is 
a systematical tool, higher-order calculations of hard kernels are demanded. In this paper we shall elucidate the 
framework for these calculations, deriving the next-to-leading-order (NLO) hard kernel for the scattering process 
TPy* ^ 7 as an example. The point is that partons in both the quark diagrams from full QCD and the effective 
diagrams for the pion wave function, carrying the momentum k — (fc + , 0, k-r), are off mass shell by k\. The difference 
between the two sets of diagrams defines the hard kernel in kx factorization theorem, a procedure similar to the 
derivation of Wilson coefficients in an effective field theory. This is the way to obtain a fcr-dependent hard kernel 
without breaking gauge invariance, since the gauge dependences cancel between the above two sets of diagrams. A 
physical quantity is expressed as a convolution of a hard kernel with model wave functions, which are determined 
by methods beyond a perturbation theory, such as lattice QCD and QCD sum rules, or extracted from experimental 
data. A gauge-invariant hard kernel then leads to gauge-invariant predictions from kx factorization theorem. 

We emphasize that the above prescription for computing a fc-r-dependent gauge-invariant hard kernel has not yet 
been fully recognized. Several NLO calculations, which include the transverse momentum dependence via on-shell 
partons carrying k = (k + , k~ ,kx), k~ ~ k T /(2k + ), have been performed in the literature [2l|, [29|, [3(J. In these 
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calculations both quark diagrams and effective diagrams are gauge-invariant, and so are hard kernels. However, 
the considered parton momentum is not a configuration described by the nonlocal matrix elements associated with 
fey-dependent hadron wave functions, because the minus component k~ should have been integrated out. Another 
subtlety is that the NLO hard kernel for the process tt^* — * 7 obtained in the above formalism turns out to be fcy- 
independent [301 ] . The parton transverse degrees of freedom in the pion wave function are then integrated out, and the 
formalism reduces to collinear factorization theorem. Moreover, we shall explain that the additional nonperturbative 
soft function introduced in [30] is not necessary for fey factorization theorem, since the infrared logarithms can be 
absorbed into the pion wave function completely. 

As stated before, the light-cone singularities [l9[ in the naive definition for fey-dependent hadron wave functions 
must be regularized. These singularities, not present in the quark diagrams, are not physical. If not regularized, 
higher-order hard kernels, computed as the difference of the quark diagrams and the effective diagrams, will be 
divergent. In this paper we shall adopt the modified definition, in which the Wilson lines involved in the nonlocal 
matrix elements for hadron wave functions are rotated away from the light cone. After the subtraction of the 
singularities, a hard kernel depends on regularization schemes unavoidably, which can, nevertheless, be regarded as 
part of the factorization-scheme dependence. This dependence, usually minimized by adhering to a fixed prescription 
for deriving hard kernels, does not cause a problem. The removal of the light-cone singularities from wave functions 
and the gauge invariance of hard kernels are the two essential ingredients for making physical predictions from fey 
factorization theorem. 

We shall demonstrate that the higher-order quark diagrams for 717* — > 7 generate two types of double logarithms, 
ln 2 (Q 2 /fcy) and In 2 x, Q 2 being the large momentum transfer squared, from the loop correction to the virtual photon 
vertex. The former does not appear in collinear factorization theorem, but the latter does [H], HH ■ It is found that 
the effective diagrams reproduce the same double logarithm ln 2 (Q 2 /fc|,), which is then absorbed into the pion wave 
function, and organized by fey resummation 0, H, H3, HH| ■ The remaining double logarithm In 2 x can be absorbed 
into the jet function, and organized by the threshold resummation [33| . Eventually, the hard kernel is free of any 
double logarithm, and its perturbative expansion is improved. It will be shown that the NLO corrections are only 
few-percent for the pion transition form factor involved in the scattering process n^f* — > 7, if setting the factorization 
scale to the momentum transfer. 

In Sec. II we calculate the 0(a s ) quark diagrams from full QCD, the 0(a s ) effective diagrams for the pion wave 
function, and the 0(a s ) jet function, and then take their difference to obtain the 0(a s ) hard kernel for nj* — * 7 
in fey factorization theorem. The gauge invariance of the fey-dependent hard kernel is proven to all orders in a s by 
induction in Sec. III. Section IV is the conclusion. 



II. 0(a s ) k T FACTORIZATION 

In this section we set up the framework for computing the hard kernel for the pion transition form factor in fey 
factorization theorem. The momentum Pi of the pion and the momentum P 2 of the out-going on-shcll photon are 
chosen as 

Px=(P+,0,0y), P 2 = (0,P 2 -,Oy). (1) 

The LO quark diagram, in which the anti-quark q carries the on-shell fractional momentum fc = (xP^~ , 0, Oy) and the 
internal quark carries P2 — fc, leads to the amplitude 

C (0) (x Q 2 ) - ir ^^ 2 ~ /Pl75 ] - t l\l ^ l75 ] (2) 

(P 2 - fc) 2 xQ 2 

with the leading spin structure Pi 75 of the pion and Q 2 = 2P\ ■ Pi. We have suppressed other constant factors, such 
as the electric charge, the color number, and the pion decay constant, which are irrelevant in the following discussion. 
The trivial factorization of Eq. ([2]) reads 0] , 



G (0) (x,Q 2 ) = J dx'd 2 k' T ^\x;x\k' T )H^{x\Q 2 ,k' T ) 
$ (0 \x;x',k' T ) = S(x - x')S(h' T ) , 



//■"V.^.A,, = - M ^r^ l7] . (3) 



Once we concentrate on the small x region, the treatment of the parton fey differs from that in collinear factorization 
theorem: fc T in the denominator of Eq. ([3]) is not small compared to xQ 2 , and the internal quark propagator should 
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not be expanded into a power series in k 2 -, [2y, [3J|. kx in the numerator, being power-suppressed by 1/Q, is combined 
with three-parton meson wave functions to form a gauge-invariant set of higher-twist contributions as in collinear 
factorization theorem. This special treatment of the parton kx characterizes the distinction between kx and collinear 
factorizations [28j |. Because of the zeroth-order wave function oc 5(k' T ), the LO hard kernel does not depend 
on the parton transverse momentum actually. 

The 0(a s ) quark diagrams corresponding to Eq. ([2|) from full QCD arc displayed in Fig.Q] in which the upper line 
represents the q quark. The factorization of the collinear divergences from these radiative corrections is referred to 

G (1) (x,Q 2 ) = J dx'd 2 k' T \^ {1 \x-x' 1 k' T )H {0 \x\Q 2 ,k' T ) + ^ \x-x\k' T )H {1 \x\Q 2 ,k' T )] [ , (4) 

where the 0(a s ) effective diagrams $W are defined by the leading-twist quark-level wave function @,[35[ 

<S>(x;x',k T ) = f ^p^e- ix ' p ^+ ik '^(0\q(y)W y (n)U n;y;0 Wo(n) A-75<?(0)kA - k)q(k)) , (5) 

with y = (0, y~,yx) being the coordinate of the anti-quark field q, ri- = (0,1, Ox) a null vector along P2, and 
\q{P\ — k)q(k)) the leading Fock state of the pion. 
The factor W y (n) with n 2 ^ denotes the Wilson line operator, 



W y {n) = Fexp 



ig / dXn ■ A(y + An) 







(6) 



The two Wilson lines W y (n) and Wo(n) are connected by a link I n ;y,o at infinity in this case 0, Equation (J5]) 
contains additional collinear divergences from the region with a loop momentum parallel to n_, as the Wilson line 
direction approaches the light cone, ie., as n — > n_ [l9(. It will be shown that n 2 serves as an infrared regulator for 
the light-cone singularities, and that the wave function depends on the additional scale £ 2 = 4(n- Pi) 2 /\n 2 \, ie., on the 
external kinematic variable. Besides, $ also depends on the factorization scale /Xf , which is not shown explicitly. Note 
that Eq. ([5|) does not reduce to the distribution amplitude in collinear factorization theorem directly, when integrated 
over kx, but a convolution of a hard kernel with the distribution amplitude [3<| . 

With one-gluon exchange, the outgoing partons from ie., the partons participating the hard scattering, carry 
the transverse momenta, so that in Eq. ((4]) depends on k' T nontrivially in the first-order factorization. Being 
convoluted with <&(°), the partons entering the NLO hard kernel ifW are still on-shell. To acquire the nontrivial kx 
dependence, must be convoluted with the higher-order wave functions $W, i > 1: the gluon exchanges in $W 
render the incoming partons of H^~\ ie., the incoming partons of the quark diagrams G^ and the effective diagrams 
$W off-shell by k T We thus derive H^(x, Q 2 , kx) according to the formula 

H^(x,Q 2 ,k T ) = G {1 \x,Q 2 ,k T )~ J dx'd 2 k T ^(x,k T -,x',k' T )H^(x',Q 2 ,k T ) , (7) 

where $W(:e, kr', x', k' T ) is defined by Eq. |(SJ) but with the q quark momentum k — (xP^ ,0,kr). As stated in the 
Introduction, the gauge dependences of G^ 1 ' and cancel in the above expression, such that H^(x, Q 2 , kr) turns 
out to be gauge- invariant. 



A. Quark Diagrams 

The loop integrals associated with the 0{a s ) quark diagrams in Figs. IHa)-(f), where the q quark carries the 
momentum k — (xP^~ , 0, kr) and the q quark carries k = Pi — k, are written, in the Feynman gauge, as 



,k T ) = —g 2 C F n 



2c 
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G^(x,Q 2 ,k T ) 



-9 C F n 



2 c 



G^(x,Q 2 ,k T ) = -ig 2 C F[ i 



2c 
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(9) 
(10) 



4 



Pi -k 



(a) 



(6) 



(c) 





(<0 L (e) t- (/) 

FIG. 1: 0(a s ) quark diagrams for it"/* — > 7 with x representing the virtual photon vertex. 



G ( j\x,Q 2 ,k T ) 



2 ,k T ) = ig 2 C F ^ 
2 ,k T ) = ig 2 C F ^ 



j4-2< 
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(11) 

(12) 
(13) 



The coefficients 1/2 in Eqs. §8§ and © arise from the definition of the self-energy corrections to external particles. 
C F is a color factor, and fi the renormalization scale. 

We work in the dimensional reduction [38j ] to simplify the calculation and to avoid the ambiguity from handling 75 
in arbitrary dimensions. The results for the self-energy corrections are 
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xQ 2 + 





H^(x,Q 2 ,k T ) , 



(14) 
(15) 
(16) 



where 1/e denotes the ultraviolet pole, and 7^ is the Euler constant. Since the external partons are off-shell by k^, 
the collinear divergences in Figs.QJa) an d 1(b) are represented by the infrared logarithms Infcf. in Eqs. (fT4"]) and (|15p . 
respectively. The internal quark in Fig.QJc) is off-shell by the invariant mass squared xQ 2 + k^, which then replaces 
the argument k\ in the infrared logarithm. 

In the small x region we drop terms suppressed by powers of x or k^/Q 2 . The loop correction to the virtual photon 
vertex gives 



G [ p{ Xl Q 2 M) = 




In 2 — 2 In — g- In ■ 



Q 2 



+2 In ■ 



In 




H^{x,Q 2 ,k T ) 



(17) 



At small x the q quark in Fig. [ljd) is energetic, implying the existence of the collinear logarithmic enhancement 
ln(Q 2 / k 2 ^) , and the internal quark is close to mass shell, implying the soft enhancement \n[Q 2 /(xQ 2 + fc 2 -)]. Their 
overlap then leads to the double logarithm \a(Q 2 /kj) ln[Q 2 /(ccQ 2 + k 2 -)] in Eq. ijTTf . In the region with x ~ O(l), 
the internal quark becomes off-shell by 0(Q 2 ), the soft enhancement disappears as ln[Q 2 /(xQ 2 + fc^)] ~ 0(1), and 

the double logarithm reduces to a single logarithm. The result of clearly exhibits the transition of the double 
logarithm in the small x region to the single logarithm in the large x region. 
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The above double logarithm deserves more discussion, which can be reexpressed as 

01 Q 2 -i Q 2 i 2 Q 2 i 2 Q 2 . i 2 ^Q 2 + &t / 1Q i 

- 2ln TT ln :E0 2 + fc 2 = ~ ln IT" 111 x0 2 +fc 2 + ln P ■ ( 18 ) 

Ihrj^ Ju\cg | f\jrj~i f\jrp Jb\c^/ \ IXrp fVj-i 

The first term is known as the Sudakov logarithm 0, [24j], which will be absorbed into the pion wave function as 
stated before. The Sudakov effect from resumming this double logarithm suppresses the contribution from the small 
kx region, ie., the region with a large impa ct p arameter The second term exists even in collinear factorization 
theorem without taking into account kx |3ll 

If, HQ 2 /(xQ 2 + k T )} ~m 2 x, which can not be factorized into the pion 
wave function. This threshold logarithm is important at small x, where the internal quark approaches mass shell. 
Hence, a jet function has been introduced to absorb In 2 x, and its resummation effect suppresses contributions from 
the small x region (33j. The third term, being of 0(1), does not require an all-order organization. 
The loop correction to the out-going on-shcll photon vertex is written as 

G^(x,Q 2 ,k T ) = ^C F Q + In 0- + In XQ2 + 4 + f ) H^(x, Q 2 , k T ) , (19) 

which does not contain a double logarithm for the following reason. In the large x region the internal quark is off-shell 
by 0(Q 2 ), and the soft enhancement disappears. In the small x region the q quark becomes soft, and the associated 
collinear enhancement is diminished by the limited phase space for the loop momentum. Therefore, there is a lack of 
overlap of the collinear and soft enhancements, and only the 0(1) single logarithm exists. 

At last, the evaluation of the box diagram Fig. QJf ) is simple, giving a power-suppressed contribution at small x. In 
the region with x ~ 0(1), ie., k + ~ 0(Q), the internal quark in Fig.QJf) is off-shell by l/[P2 '{k — I)] ~ l/Q 2 for either 
a collinear loop momentum l + ~ 0(Q) or an ultraviolet loop momentum ~ 0(Q), the same as l/(-P2 ■ k) ~ 1/Q 2 in 
the LO amplitude. Namely, the radiative correction from the box diagram does not change the LO power-law behavior, 
and its contribution is finite. In the region with small x ~ O(A), A being a hadronic scale, the LO amplitude scales 
like 1/(P 2 • k) ~ 1/(QA), while the internal quark in Fig. [TJf ) remains off-shell by 1/[P 2 • {k - I)] ~ 1/Q 2 for either 
collinear or ultraviolet I. Thus the contribution from the box diagram becomes power-suppressed and negligible, and 
we have G^\x,Q 2 ,kT) — at leading power. The above observation is consistent with the corresponding NLO 
analysis in collinear factorization theorem [3l| , which indicates the vanishing of the box-diagram contribution in the 
small x region explicitly. 

The sum of the radiative corrections from the quark diagrams Figs. []Ja)-(f) gives 

/ 

G^{ Xl Q 2 ,k T ) = Y,Gf\x,Q 2 M) 

i—a 

It is observed that all the ultraviolet poles cancel and the [i dependence disappears completely, a consequence of the 
conservation of the current that defines the pion transition form factor. It will be demonstrated in the next subsection 
that the effective diagrams for the pion wave function generate the same infrared logarithms lnfc^. 



B. Effective Diagrams 

We first explain the appearance of the nonphysical light-cone divergences in the naive definition for fc-r-dependent 
hadron wave functions. To factor out the collinear gluons in Figs. 1(d) and 1(e), the following approximation for the 
product of the two internal quark propagators has been employed [35( , 
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xQ 2 



(x-l+/P+)Q 2 + \k T -l T \ 



(21) 



where 2P£ comes from the contraction of /P 2 and 7^ in the numerators of Eqs. (| 1 1 1) and (fT2|) . The factor n^/n„ ■ I 
is exactly the Feynman rule associated with the Wilson line along the light cone, which is necessary for the gauge 
invariance of the nonlocal matrix element in the pion wave function. The first (second) term in the above splitting 
corresponds to the case without (with) the loop momentum I flowing through the hard scattering. It is easy to 
see that the right-hand side of Eq. pTjl is well-defined in the n_ ■ Z = Z + — ^ limit, if the transverse momenta k\ 





W (e) (/) (a) 

FIG. 2: 0(q s ) effective diagrams for the pion wave function. 



and |k-r — It\ 2 are dropped. That is, collinear factorization can be made gauge-invariant and free of the light-cone 
singularities. However, singularities from l + — > are developed, when the transverse momenta are included, implying 
that the factorization of collinear gluons should be performed more carefully in kx factorization theorem. This is the 
reason the naive definition is modified into Eq. © with the non-light-like vector n, which makes finite n ■ I as l + — > 0. 
The explicit expressions for the 0(a s ) effective diagrams displayed in Fig. 0(a)- (g) are written, following Eq. |5]), 

as 
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(22) 
(23) 

(24) 

(25) 

(26) 

(27) 

(28) 



where n + = (1,0, Ot) is a null vector along the pion momentum Pi, and the arguments fif and £ 2 of < & < - 1 - ) are not 
exhibited for brevity. Note that the zeroth-order wave function is given by $' ' = 8{x — x')8{kr — k' T ) here. 
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We compute the convolution of <& < - 1 - ) with the LO hard kernel in Eq. Q over the integration variables x' and 
k' T , denoted by (g> below: 



(i) 



H {0) 



dx'd 2 k' T $f\x, k T ; x', k' T )H {0) (x, Q 2 ,k' T 



(29) 



The self-energy corrections in Figs.[2ja) and^b) are similar to the quark diagrams in Figs.QJa) andQJb), respectively, 
and the results are 



(i) 



8tt 



4> 



87T \ £ 



In 



In 



4^ f 2 



2)H^{x,Q\k T ) 



2) H^(x,Q 2 ,k T ) 



(30) 
(31) 



Similarly, the contribution from the box diagram Fig. [2jc) is power-suppressed in the small x region, and we have 

When evaluating Eqs. ([23]) - (f2"8")) . the sign of the plus component n + of the vector n is arbitrary, which could be 
positive or negative (n~ has a positive sign, the same as of Pj - )- Choosing n + < 0, ie., n 2 < as in (El. Isl [Ill fl2| . 
Fig. HJd) leads, in the small x region, to 
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(32) 



which reproduces the Sudakov logarithm ln 2 (Q 2 /fcy) from Fig. QJd) in Eq. (fl"8"|) . noticing the scale £ 2 = n^/n+IQ 2 . 
The light-cone divergences are regularized in the price that the universality of a wave function is lost, for it depends 
on the external kinematic variable through £ 2 . This problem can be alleviated b y ex tracting the evolution in £ 2 from 
Eq. §5§ [l9j], ie., by resumming ln 2 (C 2 //4) in Eq. ^ into the Sudakov factor 0,0. The initial condition of the 
evolution is universal, like a distribution amplitude in collinear factorization theorem. We stress that the Sudakov 
resummation, accurate up to fixed loops, does not remove the ( 2 dependence of a wave function completely. That is, 
nonfactorizability may occur at subleading level in kx factorization of the pion transition form factor. 

The hard kernel associated with $i , ie., the second term in Eq. (f2Tj) . demands the physical range of l + to be 
— k + < l + < fc + , which corresponds to the range of the parton momentum fraction 1 > x' > 0. As computing the 
convolution of with H^°\ this fact should be taken into account. Moreover, we assume ( 2 ~ Q 2 by choosing 
\n + \ ~ n~ to avoid creating the additional large logarithm ln(£ 2 /<2 2 ). The leading-power expression for Fig. &,e) is 
then given, in the small x region, by 



Air 



2 + fc T) jff( o )(a 



(33) 



where terms vanishing with k 2 ^ — > have been dropped. It is found that Fig. [2lje) does not generate a large double 
logarithm with £ 2 ~ Q 2 . 

It is interesting to obtain the results corresponding to n + > for Figs. [2^d) and[2][e). One simply analytically 
continues Eqs. (|32[) and (|33[) into the region with n 2 > by means of the principle-value prescription, 
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We then derive from Eq. 
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(35) 



which can be confirmed by calculating the loop integral in Eq. (|25ll directly for n 2 > 0. It shows that the choices 
n 2 > and rt 2 < lead to expressions different only by a constant term. Because the n-dependent double logarithms 
cancel in the summation 
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(a) (b) 
FIG. 3: 0(a s ) diagrams for the jet function. 



Applying the variable change I — > — Z, and the transformation n — > — n and k — > fc, Eq. (|27[) becomes identical to 
Eq. (|25[) . Therefore, the result from Fig. [5Jf) is the same as of Fig. HJd), but with the replacement of k rs Pi by fc, 
ie., £ by x£. Keeping terms which do not vanish with k 2 -, — > 0, we have 



if'- '' does not depend on the sign of n 2 actually. 



■r 4?r V e k 2 e~iv 



In 



x 2 C 2 



In 



x 2 C 2 



(37) 



where the double logarithm, being large in the region of x ~ 0(1), attenuates with the decrease of x. It should 
disappear, after combined with the contribution from Fig. [2jg) , since such a double logarithm is absent in the 
corresponding quark diagram Fig. [Ije) in any region of x. The same variable transformation relating to is 
not applicable to $ 9 1 ' > , for the latter involves the nontrivial convolution with H^°\ Hence, <f>^ <S> is expected 

to have an expression different from ® H^°\ Retaining terms which are finite as kx — > 0, Fig.[^g) leads, in the 
small x region with xQ 2 3> x 2 ( 2 , to 



47T 



CWn 



H^(x,Q 2 ,kr) 



(38) 



The cancellation of the double logarithms in the summation of Eqs. (|37|) and |38|) is obvious. For a similar reason, 



+ $g ) (8 -ff^ is independent of the sign of n 2 . 
Summing all the above 0{a s ) quark-level wave functions, we derive 



<&« ® # (°) = g$ 



(i) 



47T 



In 



, 2 ? 



2/-2 



In 



In 



x\_ 



2tt 5 



21,2 



Q 2 k 

H^{x,Q 2 ,k T ) 



(39) 



In contrast to Eq. (|20p . which is independent of the renormalization scale n, the above expression depends on the 
factorizations scale /if. The Sudakov resummation and the renormalization-group method can be applied to organize 
the logarithms ln 2 (£ 2 /fcf,) and ln(^ 2 /fc^) to all orders, respectively @]. 



C. 0{a s ) Hard Kernel 



We renormalize Eq. (|3"9"|) in the modified minimal subtraction scheme, and then take the difference of Eqs. (piO)) 
and ([39)1 to obtain the 0(a s ) hard kernel for the pion transition form factor. It is easy to find that the hard kernels 

H<-°\ and 



o-U) - r (.l) 



(i) 



G 



(i) 



y = v_rj — yy u = associated with Figs. QJa)-(f) are all free of the infrared logarithms lnfef, as claimed 
before. Compared to (3(j, we do not need the additional soft function S to achieve this cancellation. The difference is 
that the self-ener gy c orrections to the Wilson lines have been included into the set of effective diagrams for the pion 
wave function in [301 ] . Hence, S must be introduced to remove these artificially included infrared divergences. We 
stress that the self-energy corrections to the Wilson lines do not exist, because such diagrams are not generated in 
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the derivation of factorization theorem using the diagrammatic approach 0] • This observation is consistent with the 
postulation that the gauge fields appearing in the Wilson lines in Eq. © are regarded as bare fields [l9[ . 

After subtracting the effective diagrams from the quark diagrams, the resultant hard kernel depends on the fac- 
torization scheme that defines the renormalization of Eq. (|3D[) . The quark diagrams do not have such a scheme 
dependence as shown in Eq. (|20[) . When making a physical prediction from factorization theorem, one convolutes the 
hard kernel with a model for the pion wave function (not with the effective diagrams), so that the scheme dependence 
in the hard kernel remains. As stated in the Introduction, the scheme dependence of physical predictions is usually 
minimized by adhering to a fixed prescription for deriving hard kernels, which will be elucidated below. The sum of 
the 0(a s ) hard kernels is written as 

/ 

H^(x,Q 2 ,k T ) = ]T^ (1) (:r,Q 2 ,fcT) 



i— a 



I 2 /-2 

■j—Cp I -In * l , 2 +21n)-ln — ^ ;2 -In' " , 2 
47r \ xQ 2 + k\ Q 2 xQ 2 + fcf xQ 2 + k\ 



(40) 



The Sudakov logarithm ln 2 (Q 2 /fc§,) in Eq. (JT5J! for has been cancelled by that in Eq. §^ for <^ 1} ® H^°\ but 
the threshold logarithm In 2 [Q 2 / (xQ 2 + fcf,)] remains in H^ l \ The large threshold logarithm can be absorbed into a 
jet function [33[, so that the pertubative expansion of the hard kernel is further improved. At small x, a collinear 
enhancement arises from the region with the loop momentum parallel to the internal quark momentum P 2 — k ~ P-i. 
To factorize this collinear gluon into the jet function, we replace the q quark by the eikonal line in some direction u 
[40l |4]| as shown in Fig. [3Ja). Similarly, we choose u 2 ^ to avoid other infrared divergences, such as those from 
I parallel to Pi, which have been absorbed into the pion wave function. Including the self-energy correction to the 
internal quark, Fig.[3](b), we arrive at the complete set of diagrams for the jet function at 0(a s ). 

Figure [3] has been evaluated in [33| . focusing only on the double- logarithm piece In 2 x. Here we work out the 

single-logarithm and constant pieces too. The explicit expression of the loop integral J<j associated with Fig. [3]Ja) is 
referred to [33J. We obtain, for u 2 < 0, 

tf >„<«> = gov (1 + 1„ ^ - 1" 2 ^ + m + > - £ ) flP. M , («> 

with the scale ( 2 = 4(u • P 2 ) 2 /|w 2 |. Figure [31(b) gives a result identical to Eq. (Til)]) for Fig. QJc): 

Note that the sum J (1) = J^ 1} + J b (1) is free of ultraviolet divergences and /x-independent. That is, the factorization of 
the jet function does not modify the renormalization-group behavior of the hard kernel. As expected, the jet function 
is characterized by the invariant mass of the internal quark. 

Define £ 2 = vQ 2 and £ 2 = v u Q 2 with v and v u being constants of O(l). The hard kernel, after subtracting the 
0{a s ) jet function, is given by 

[H/J]^(x,Q 2 ,k T ) = H^(x,Q 2 ,k T ) - J^(x,Q 2 ,k T )H^(x,Q 2 ,k T ) 



ln r n 2 + k 2 ~ 2(ln V + ln Vu) ln rO 2 + k 2 + 2 ^ 



3 — In v u + In v u + 2 ln ^ 



H^(x,Q 2 ,k T ), (43) 



in which the double logarithms have been completely removed. Different values of v and u u correspond to different 
factorization schemes. Adopting v — 1, ie., ( 2 = Q 2 as in @, and v u — \, Eq. (|4"3"1) reduces to 



[if/J]( 1 )(a;,Q 2 ,fc T )=-? 1 CF In , 2 +21nx + 3-^ gW(x,Q 2 ,fc T ). (44) 
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(6) 

FIG. 4: (a) Diagrams for XdG N+ ^ /dX, where the bubbles represent , and the squares contain the special vertex v a . (b) 
Diagrams for Xd&^/dX. 



Employing the factorization scale m = Q and the asymptotic model of the pion wave function, the same as in the LO 
analysis in fcy factorization theorem 42], the NLO corrections are found to be only 5%. That is, the NLO corrections 
are not expected to affect much the LO results for 7T7* — > 7. Our conclusion is drawn under the specific factorization 
scheme with v = u u = 1. It requires an examination whether NLO corrections are also negligible under the same 
scheme in other exclusive processes containing pions, such as the pion form factor involved in 7r7* — » tt. 



III. GAUGE INVARIANCE 



In this section we prove the gauge invariance of the fey-dependent hard kernel for the pion transition form factor 
by induction. We first show that the kr factorization constructed in the Feynman gauge @, [H[ holds in an arbitrary 
covariant gauge d ■ A — with the gauge parameter A, in which the gluon propagator is given by (~i/l 2 )N^ u (l) with 
the tensor 



IH V 
1 

It has been argued that the replacement 



N^(l)=g^-(1-X)— . (45) 



n^F n I" 

9^^— ,=9^^, (46) 

for a collinear gluon propagator in the Feynman gauge extracts collinear divergences correctly @, [3|| . In the arbitrary 
covariant gauge we just need to modify the above replacement into 

N^(l) - — (1 - A) Li- = N^(l)^L , (47) 

and then the procedures for deriving factorization theorem in 0, [35| follow: the Ward identity is applied to all the 
contractions of l v , leading to the factorization of the collinear gluon. The factor n_njn_ ■ I explains how the Wilson 
lines are generated in factorizing hadron wave functions. For more details, refer to [71 [35|. 
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The kr dependence in a hard kernel implies that the partons entering the quark diagrams and the effective diagrams 
for the pion wave function are off-shell by kj.. The LO hard kernel H^°'{x, Q 2 , kr) in Eq. ([3]), which does not contain 
a gluon, is independent of the gauge parameter A. Beyond LO, the gauge invariance of a hard kernel is a consequence 
of the gauge-dependence cancellation between the above two sets of diagrams. Assuming that the hard kernels defined 
by 



H^{x,Q 2 ,k T ) = G u) (x,Q 2 ,k T )-J2 I dx'd 2 k' T ^{x,k T ;x',k' T )H^- i \x',Q' 1 ,k' T ) , (48) 

i— 1 

are gauge-invariant for j = 1, 2, • • • JV, we shall prove the gauge invariance of the 0(a^ +1 ) hard kernel 

N+l . 

H^ N+1 \x,Q 2 ,k T ) = G^ N+1 Hx,Q 2 ,k T )-y / dx l d 2 k l T ^{x,k T \x',k l T )H^ N+1 ^\x',Q 2 ,k' T ) , (49) 



i=l 



using the method proposed in [43j |. Note that the external quark spinors in the nonlocal matrix element in Eq. ([5|) 
have absorbed half of the self-energy corrections. Another half goes into the higher-order wave functions, giving the 
coefficients 1/2 in Eqs. (f2"2")l and (|2"5|) . The same explanation applies to the appearance of 1/2 in Eqs. ([5]) and ^ for 
the 0{a s ) quark diagrams. To discuss the gauge dependence, we consider the full self-energy corrections to the quark 
diagrams G and to the effective diagrams $. 

Applying the differential operator \d/d\ to H^ N+1 \ it acts only on the gluon propagators in G^ N+1 ^ and $W on 
the right-hand side of Eq. (|4*9"| . leading to 

X-^-N^ = X 1 -^- =vJl^N a " + N^ a l w ) , (50) 
dX l z 

with the special vertex v a — l a /(2l 2 ). The derivatives XdH^ N+1 ~^ /dX vanish due to the gauge-invariant assumption 
associated with Eq. (|48|) . The loop momentum Z M (Z") in Eq. (|50| contracts with vertices in the diagrams of G^ N+1 ^ 
and $W, which are then replaced by the special vertex v a . Summing all the quark diagrams with various differentiated 
gluons and employing the Ward identity, only those, in which the special vertex is located at the outer ends of the 
valence quark lines, are left [35j,|43( as shown in Fig.[3Ja). These diagrams come from the second terms in the following 
Ward identities associated with the quark and the anti-quark, respectively, 

A75H = ^75+ At* > (51) 

where /Pi 75 is the leading spin structure appearing in the expressions for the quark diagrams. We have the similar 
Ward identities for the effective diagrams with /1+75 being substituted for /Pi 75. The first term is cancelled by one 
of the two terms from the contraction of I with the adjacent vertex. If all the external quarks were on mass shell, 
ie., kr = 0, the second terms also vanish due to jk /Pi = /Pi ]k = 0, implying XdG^ N+1 ^/dX = and Xd^^/dX = 0. 
That is, the quark diagrams from full QCD and the effective diagrams for the wave function with on-shell partons are 
gauge- invariant . 

For the differentiated quark diagrams G^ N+1 ' in Fig. HJa), the gluon emitting from the special vertex attaches all 
the lines inside G^ N \ Adopting Eq. (|47|) and the procedures in [35|], XdG^ N+1 ^ / dX is factorized into the convolution 
of (JW with the differentiated <&( 1 ) at leading power in 1/Q: 



X^-G^ N+l \x,Q 2 ,kr)= j dx'd 2 k' T X^-^{x,k T ;x\k' T )G^ N \x',Q 2 ,k' 7 
dX J dX 



) ■ (52) 



For illustration, we display the effective diagrams for XdQ^ /dX in Fig. 2Kb) explicitly. We repeat the above steps for 
the differentiated wave function $( l \x, kx] x' , k' T ), and obtain 

x4r$ il) (x,k Tl x",k'+) = / dx'd 2 k' T x4r^ (1) {x,k T ;x',k , T )^ l - 1 Hx\k' T ;x",k^) . (53) 
dX J dX 

Combining Eqs. (fB"2")) and ([5^1) . the differentiation of the 0(a^ +1 ) hard kernel gives 

\±H( N+1 \x,Q 2 ,kT) = [ dx'd 2 k' T X^-^ 1 \x,k T ;x',k' T ) 
dX J dX 

r N r 

x GW(x',Q 2 ,k' T )-Y / dx"d 2 k'^(x',k' T ;x" > k'^)H( N - i \x",Q 2 ,k^) , (54) 



2 = 
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where the term in the square brackets diminishes because of Eq. (|48f for j = N. We then prove the gauge invariance 
of the 0(a^ +1 ) hard kernel. The hard kernels and the resultant predictions from the kx factorization theorem are 
thus gauge-invariant to all orders by induction. 

IV. CONCLUSION 

In this paper we have elucidated the framework for the higher-order calculations in kj~ factorization theorem, which 
is appropriate for QCD processes dominated by contributions from small momentum fractions. The point is that 
partons in both the quark diagrams from full QCD and the effective diagrams for hadron wave functions are off mass 
shell by k\. Their difference gives the gauge-invariant /cT-dependent hard kernels, since the gauge dependences cancel 
between the two sets of diagrams. The gauge invariance of the hard kernels for the scattering process 7r7* — > 7 
in kx factorization theorem has been proven to all orders by induction. The proof can be easily generalized to 
other processes. We have explained that the light-cone divergences in a naive definition of fc^-dependent hadron 
wave functions are regularized by rotating the Wilson lines away from the light cone. This procedure introduces a 
regularization-scheme dependence, which, however, can be regarded as part of the factorization-scheme dependence, 
and minimized by adhering to a fixed prescription for deriving hard kernels. The gauge invariance of a hard kernel 
and the removal of the light-cone singularities are the two essential ingredients for making physical predictions from 
kr factorization theorem. 

We have calculated the NLO ^-dependent hard kernel for 7r7* — > 7 in the region with a large momentum transfer 
Q 2 and a small momentum fraction x. We have demonstrated that the infrared logarithms \a.k T , reflecting the 
collinear divergences, cancel between the quark diagrams and the effective diagrams exactly. Hence, there is no 
need to introduce the additional nonperturbative soft function in kx factorization theorem. The quark diagrams 
generate the double logarithms ln 2 (Q 2 /k T ) and In 2 a; from the loop correction to the virtual photon vertex. It has 
been shown that the former is absorbed into the pion wave function, and the latter into the jet function, confirming 
the observations made in our previous works [1, [33j]. Note that the factorization of the jet function does not alter 
the renormalization-group behavior of the hard kernel. Eventually, the NLO corrections in ttj* — > 7 amount only to 
5% under a specific factorization scheme with the factorization scale set to the momentum transfer. NLO corrections 
under the same factorization scheme in other processes, such as the pion form factor and heavy-to-light transition 
form factors, will be examined elsewhere. At last, we mention that the off-light-cone effects from yr 7^ have been 
found to be sizable in some heavy-to- light correlators based on a QCD-sum-rule analysis recently [44(. 
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of Physics, Academia Sinica and Department of Physics, National Taiwan University for their hospitality during his 
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